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Abstract 

The young star clusters we observe today are the building blocks of a new generation of 
stars and planets in our Galaxy and beyond. Despite their fundamental role we still lack 
knowledge about the conditions under which star clusters form and the impact of these 
often harsh environments on the evolution of their stellar and substellar members. 
We demonstrate the vital role numerical simulations play to uncover both key issues. Using 
dynamical models of different star cluster environments we show the variety of effects stellar 
interactions potentially have. Moreover, our significantly improved measure of mass segre- 
gation reveals that it can occur rapidly even for star clusters without substructure. This 
finding is a critical step to resolve the controversial debate on mass segregation in young 
star clusters and provides strong constraints on their initial conditions. 



1 Introduction 



According to current knowledge, planetary systems form from the accretion discs around 
young stars. These young stars are in most cases not isolated, but are part of a cluster 
(e.g. Lada & Lada 2003; Evans et al. 2009). Densities in these cluster environments vary 

(e.g. Tj Chameleontis) to 10^ pc~^ (e.g. Arches 



considerably, spanning a range of 10 pc 
Cluster). Though it is known that discs disperse on a time-scale of 1-10 Myr (Haisch et al. 
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200 H |Hillenbrandl|2002t |Sicilia-Aguilar et al.| [20061 [Currie et aL]|2008[ ) and t hat in dense 
clusters {n > 10"^ pc"'') the disc frequency seems to be lower in the core (e.g. Balog et al. 



2007), it is an open question as to how far interactions with the surrounding stars influence 
the planet formation in clusters of different densities. An encounter between a circumstellar 
disc and a nearby passing star can lead to a significant loss of mass and angular momentum 



from the disc. While such isolated encounters have been studied in a large variety (Heller 



1993 


Clarke & Pringle||1993[ Ostriker||1994| Heller||1995| Hall et al.||1996 Han||1997 


Pfalzner 


2004 


Pfalzner et al.|2005t 


Moeckel & Bally|2006 20071 


Kley et al.|2008 


) , only a few numerical 



studies have directly investigated the effect of stellar encounters on circumstellar discs in a 



dense cluster environment (Scally & Clarke 2001; Adams et al. 2006). 



The dynamical evolution of a star cluster leaves a variety of imprints in the phase space 
of its stellar population which are good tracers of the dynamical age of the cluster. One of 
the most widely discussed aspects is that of mass segregation. Due to energy equipartition 
- hence via two-body encounters - the more massive particles tend to settle towards the 
cluster centre over time while the lower-mass particles are preferentially pushed to the outer 



parts (Spitzer 1969, Farouki et al. 1983 Spurzem & Takahashi 1995; Khahsi et al. 2007). 



However, it is a much more challenging task to identify mass segregation observationally in 
real objects than theoretically from 'clean' numerical simulations. The investigation of mass 
segregation in young stellar systems is of particular interest for a deeper understanding of 
the star formation process. 



2 Star-Disc Encounters in Young Star Clusters 



2.1 Method 



We follow the idea of Scally Sz Clarke (2001), combining a simulation of the dynamics of 



a cluster to determine the interaction parameters of close encounters between stars in the 
cluster with results from studies of isolated star-disc encounter simulations. Throughout 
this work we assume that initially all stars are surrounded by protoplanetary discs. This is 
justified by observations that reveal disc fractions of nearly 100 % in very young star clusters 
(e.g. [Haisch et aL]|2000t |Lada et aLpOOOj [Haisch et aL]|2001[ |Hillenbrand]|2005D . 

The dynamical cluster models contain only single stellar components without consider- 
ing embedded gas. The simulations were performed with nbody6++, nbody6-GPU (Spurzem 



1999 



2003 



Aarseth 2003 2008), and starlab (McMillan 1996 Portegies Zwart et al. 2001 Hut 



. The details of the numerics can be found in (Olczak et al. 2006, 2011, in prep). 
The encounter-induced transport of mass and angular momentum in protoplanetary 



discs is calculated via fit formulae based on simulations of star-disc encounters (Olczak et al. 



2006 Pfalzner et al. 2006 Pfalzner Sz Olczak 2007). These quantify upper limits due to the 



restriction to low-mass discs and co-planar, prograde, and parabolic orbits, which are the 
most perturbing. A simplified prescription terms stars that have lost more than 90 % of their 
initial disc mass as "discless" . The results presented here apply to a scenario where all stars 
have initially a disc size of 150 AU. A lerger range of models has been investigated in the 
cited publications. 
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Figure 1: Left: Fraction of star-disc system over time in the entire ONC {R = 2.5 pc) and 
its dense core, the Trapezium cluster {R = 0.3 pc). Right: Average cumulative number of 
encounters vs. stellar mass after 1 Myr of dynamical evolution. 



2.2 Results 

From a numerical model of the Orion Nebula Cluster (ONC) we demonstrate that the 
encounter-induced disc-mass loss becomes significant in its dense core, known as the "Trapez- 
ium Cluster" (Fig. [T] left). Up to 20% of the discs can be destroyed by gravitational inter- 
actions within a radius of 0.3 pc (compared to ~5% in the entire cluster of 2.5 pc radius). 

The massive stars in the centre of such a stellar cluster act as gravitational foci for 



the lower mass stars (Pfalzner et al. 2006). This becomes evident from the average number 



of encounters as a function of stellar mass (Fig. [T| right): the number of disc-perturbing 
interactions is nearly constant for low- and intermediate-mass stars but increases largely for 
high-mass stars. Because discs are most affected when the masses of the interacting stars are 



unequal (Olczak et al. 2006 Moeckel & Bally 2007) massive stars dominate the encounter 



induced disc-mass loss in star clusters like the ONC. For low-mass stars the mass- loss occurs 
through few strong encounter events, whereas the disc of high-mass stars is removed via a 
steady nibbling by many encounters with stars of lower mass. 

The dynamically outstanding role of high-mass stars does also affect the stellar multi- 
plicity in the ONC. It turns out that the most massive star (M* = 50 M©) has on average 
of the order of 200 so-called "capturing encounters" - encounter events with an eccentricity 
e < 1 - in the first 5 Myr while a star with = 1 has less than one. These interactions 
lead to the formation of relatively stable configurations termed as "transient bound systems" 
(TBS). The left plot in Fig. [2] shows a striking increase of tbound towards more massive stars 
by more than one order of magnitude. These TBS translate directly into an observed appar- 
ent binarity (Fig. [2| right): of the five most massive stars in the ONC on average 30 % would 
appear the be in a binary system at the cluster's current age of ~1 Myr. The average for all 
13 OB stars in the ONC would be more than 10 %. 

The encounter-induced disc-mass loss becomes much more apparent in more extreme 
environments. Using a numerical model of the Arches cluster, one of the densest and most 



massive clusters in our Milky Way (Figer et al. 



current age of ~2.5Myr (Najarro et al. 2004; Martins et al. 2008) nearly all circumstellar 



2002 



Stolte et al. 2005), we find that at its 
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simulation time [IVIyr] 



Figure 2: Left: Average duration of a bound state, tbound, as a function of the primary mass 
Ml. Right: Fraction of the 5 (red) and 13 (green) most massive stars that are part of an 
encounter-induced binary as a function of cluster age. The lines are Bezier curves of the data. 




Figure 3: Left: Fraction of star-disc system over time in the entire Arches cluster. Right: 
Average disc fraction vs. stehar mass after 2Myr of dynamical evolution. 



discs could have been destroyed just due to gravitational interactions (Fig. [3| left). A striking 
result is that the low fraction of members that could retain their discs to some degree at an 
age of 2 Myr is populated by stars with masses between 2M0 and 30 (Fig. |3) right). 



3 Mass Segregation in Young Star Clusters 



3.1 Method 



As a proxy for mass segregation we extend the method A^mst developed by Allison et al. 



p009l ) (see also [Cartwright fc Whitworth|[2004| [Schmeja fc Klessen|[2006l ) . In summary, the 
authors use the minimum spanning tree (MST), the graph which connects all vertices within 



a given sample with the lowest possible sum of edges and no closed loops (Gower & Ross 
1969). The length of the MST, /msT; is a measure of the concentration or compactness of 
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Figure 4: Left: Tmst and Amst for the 5, 10, 20, 50, 100, 200, 500, 1000 most massive 
stars of an initially mass segregated cluster of Ik stars and S = 0.3. The error bars and 
line thickness mark the la, 2a, and 3a uncertainties. Right: Fmst of the 5, 10, 20, 50, and 
500 most massive stars over time for a cluster with Ik members and cold initial conditions 
{Q = 0.1). The filled regions indicate la uncertainties. 



a given sample of vertices. Mass segregation of a stellar system of size N is quantified by 
comparing /mst of the n most massive stars, ^msT' with the average Zmst of k sets of n 
random cluster stars, < /jlgx >i and its standard deviation, AZ^Igrp. 

Our method A^mst involves a crucial modification of A^^grp that boosts its sensitivity: 
we do not use directly the sum of the edges /mst as a measure yet their geometric mean 7msTj 



/ n s l/n 

7MST = ( n ^0 = exp 



1 

— > In Cj 

1=1 



(1) 



and its associated geometric standard deviation A7mst5 where ei are the n MST edges. We 
obtain the new measure Tmst via a proper normalisation: 

ref 

^MST = -Mfi ' ^^MST = MwS,-! ■ (2) 
/MST 

The geometrical mean acts as an intermediate-pass that damps contributions from extreme 
edge lengths very effectively, hence significantly reduces the contribution of any "outlier" . 



3.2 Results 



Subr et al. 



(2008) we have created initially mass-segregated star clus- 



Using the method of 

ters to verify our algorithm. In the left plot of Fig. |4] we compare our improved measure 
rMST ("geometric") with Amst ("arithmetic"). We find that our measure of mass segrega- 
tion, FmsT) detects an intermediate degree of mass segregation, S = 0.3, with at least 3a 
significance while Amst provides only a very weak la significance. 

We have used this highly improved sensitivity to follow the dynamical mass segregation 
of a numerical star cluster model. The initial configuration is based on our ONC model as 



described in previous publications (e.g. Olczak et al. 2010). Note that we use a spherically 
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symmetric model with a smooth density distribution without any substructure. However, 
our stellar system is initially collapsing, starting from a virial ratio Q = 0.1. The simulations 
were carried out until a physical age of 5Myr, corresponding to 13.5 N-body time units. On 
the right-hand side of Fig. [4] we see that Tmst increases rapidly within the first 2 Myr (or 
5.5 N-body time units) and then saturates. The different mass groups are clearly separated, 
inversely correlated with the sample size. 

4 Discussion and Conclusions 

From numerical simulations of star cluster dynamics and star-disc encounters we find that 
pure gravitational interactions of disc-surrounded stars can lead to a significant depletion 
of their circumstellar matter. The disc-mass loss increases with cluster density but remains 
unaffected by the size of the stellar population. In moderately dense clusters like the ONC 
it is the massive stars that dominate the encounter-induced disc-mass loss by gravitational 
focusing of low-mass stars. 

Massive stars have a much higher probability to be involved in a capturing encounter in 
such a cluster environment than solar-mass stars. In the ONC at least 10 % of the OB stars 
form a (transient) bound state due to dynamical capture that would be observed as a binary. 
The properties of the these systems strongly depend on the cluster age with a development 
towards smaller periastra, higher mass ratios and longer duration of the bound states. 

In extreme environments like the Arches cluster the pairwise interaction of all stars leads 
nearly to a complete destruction of the entire disc population. The preferential survivability of 
discs around B-type stars in the Arches cluster implies that these could be the best candidates 
for detecting planetary systems in Galactic starburst clusters. 

We have developed a significantly improved method to measure mass segregation based 
on the minimum spanning tree (MST). In general, using only the ten to twenty most massive 
stars Fmst provides a robust and sensitive measure of mass segregation for the entire popu- 
lation of star clusters in our Galaxy. In particular, very low degrees of mass segregation can 
be detected in massive clusters like NGC 3603 that consist of 10k or more stars. 

We note that mass segregation in a collapsing, intermediate-size stellar cluster of Ik 
stars can occur very quickly, i.e. within only a few crossing times. This finding demonstrates 
that rapid mass segregation (in terms of dynamical time scale) does not require substructure. 
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